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In this talk I first present a short review of fluctuations in e+e^-annihilations. I then describe some new 
results on FD correlations. Experimental analyses of pp and AA correlations indicate a very small production 
radius. This result relies very strongly on comparisons with MC simulations. A study of the approximations 
and uncertainties is these simulations imply that it is premature to draw such a conclusion from the data. 
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1. FLUCTUATIONS 

In bremsstrahlung in QED, the emission of a photon does 
not change the current for subsequent emissions. This implies 
that the photon multiplicity is described by a Poisson distri- 
bution. As the gluons carry colour charge, the emission of a 
gluon in QCD changes the current relevant for the subsequent 
emissions. An initial hard gluon will radiate many softer glu- 
ons. This leads to a cascade with an exponential growth in 
multiplicity and to large fluctuations. The multiplicity distri- 
bution is approximately satisfying KNO scaling, with a width 
which is proportional to the average multiplicity. This short 
review is divided in three parts: 

A. Multiplicity distribution for partons 

B. Effects of hadronization 

C. Fractal structures 



A. Multiplicity distribution for partons 



In e+e^-ann. the emission of a gluon from a qq pair is 
a coherent emission from e.g. a red and an anti-red colour 
charge. The result corresponds to a colour dipole with the 
following distribution (in the leading log approximation): 
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The phase space is a triangular region in (y, In A:^) -space given 
by \y\ < \n{s/k^)/2, and within this triangle the density is just 
given by a. The emission of a second (softer) gluon is given 
by two dipoles, one between the quark and the first gluon, and 
one between this gluon and the antiquark. The phase space is 
enlarged compared to that for the first gluon, and coiTesponds 
to the folded surface shown in fig. [1^. For subsequent gluons 
the phase space is further increased, and corresponds to the 
fractal surface in fig.[T]5. 

We let P{N,L — Ins) denote the distribution in the number, 
A^, of dipoles (which neglecting qq pair creation is equal to the 
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FIG. 1: The phase space for gluon emission in e+e^-ann. is a tri- 
angular region in the (>■, K = In fc^) -plane. The height of the triangle 
is given by L = Ini-. When one gluon is emitted at {yi , Ki ) a second 
gluon with smaller can be emitted in the larger phase space shown 
in fig. (a). Further softer emissions give the fractal phase space in 
fig. (b). 



number of gluons plus one). The Laplace transform fp (y, L) is 
given by 

¥{y,L)^Y.e-y''P{N,L). (2) 

N 

It is then relatively easy to derive the relation (see e.g. ref. fP]) 
d^ln^P 
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Expanding eq. (|2|i in a Taylor series in y gives 
£P ^l-y{N) + ^'f{N^) + ... 

Inserted in eq. (O this also implies the following equations for 
the average and the variance of the multiplicity distribution: 
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The solutions to these equations are Bessel functions, and for 
high energies we get 



{N) - L'/'^exp(2yotoZ) 
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We see that the width of the distribution is proportional to 
(A^), in agreement with KNO scaling, and thus for large the 
distribution is much wider than a Poisson distribution. 

The anomalous dimension is given by a logarithmic deriva- 
tive. Note that in the literature one can find different defini- 
tions, where the derivative is taken with respect either to Ini 
or to InW: 
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I want here to add a few comments: 
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• (A^) depends on the resolution, i.e. on the cut-off for soft 
emissions. It is here essential to have a cut-off which is 
defined locally, and not fixed in e.g. the overall cms. In 
the dipole cascade model this is chosen as kj^ measured 
in the rest frame of the emitting dipole. 

• The running of is important. With a constant a, the 
multiplicity would grow proportional to exp(\/dc i). 
i.e. much faster than the result in eq. (|5]). 

• Non-leading corrections are very large. Terms sup- 
pressed by factors 1 / ^/Z in the evolution equation give 
extra powers of L as factors in the solution. Also in- 
cluding NLL terms, the result is sensitive to effects of 
still higher order. For a more detailed discussion see 
e.g. ref. lH. 



B. Hadronization effects 

In string fragmentation the average hadron multiplicity, (n) , 
in a single qq system is proportional to ln(,s-/io), where sq 
is a scale of order IGeV. The variance of the fluctuations, 
V = (n^) — (n)^, is proportional to («), and the distribution 
is thus relatively more narrow for higher energies and larger 
(«) . For a system of a qq pair and a number of gluons we get a 
set of string pieces with (squared) masses 1 — {qi + qi+i)'^, 
where qi is the momentum of parton / and the partons are or- 
dered in colour. For the average multiplicity we then get 



Y^\n{si,i+i/so) 
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which just corresponds to the length of the baseline of the 
fractal surface in fig.[T]3. 

In eq. ( fTOl i the mass of a string piece should not be smaller 
than ^/sQ, as it otherwise would give a negative contribution. It 
is, however, also possible to define an infrared stable measure, 
called the ^--measure, which is insensitive to the cut-off for the 
perturbative cascade Simulations show that a high energy 
even with little gluon radiation and an event at lower energy 
but more radiation give equally many hadrons if they have the 
same A,-value. Thus the multiplicity depends only on the mea- 
sure X, which corresponds to an "effective string length". The 
probability density P{X,L) satisfies the same evolution equa- 
tion (|3]l as P{N,L), only with different boundary conditions. 



This implies that the distribution in X has the same high energy 
behaviour as the distribution in dipole multiplicity, eq. (|5]l. 

The fluctuations in the hadron multiplicity get contributions 
from both the perturbative cascade and the soft hadronization, 
and we get approximately 



V{n,s) _ V{X,s) V{n,{X)) 



(11) 



Here the first term corresponds to the fluctuations in the cas- 
cade and the second one to those in the hadronization pro- 
cess. At high energy the first contribution dominates, while 
the second is larger at energies below 40-50 GeV. In ref. IH 
it is shown that the total result satisfies KNO scaling with 
V = const ■ (n)^, although this is not the case for any of the 
contributions separately below the top LEP energy. At higher 
energies the hadronization contribution can be neglected. 



C. Fractal structures and intermittency 

The baseline of the surface in fig. [TJ) looks like a Koch 
snowflake curve. Looking at the curve with a coarser res- 
olution means that only emissions with k± > A:xres are in- 
cluded. This corresponds to cutting the surface in fig. [TJ) at 
a higher level, corresponding to Infe^^.^^. The length of this 
curve satisfies again the same evolution equation. The bound- 
ary conditions are different, and with the notations L = Ini 
and K = InA;^ the length will be proportional to the expression 

Length-^ L^^'^KrisSy^p{2^/a^-2y/aoK^s) (12) 

It is also possible to define a (multi)fractal dimension for this 
curve 



^ ^ dlnlength ^ 

c/ln(Kres) 
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(It is called a multifractal as the dimension varies with the 
resolution K,es-) 

Fluctuations in small phase space regions, e.g. slices in 
rapidity Ay, can be described by the scaled moments 
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If the total phase space is divided in M equal pieces (which 
gives Ay ~ Y /M) a scaling behaviour Fq °= M*)"'' is interpreted 
as an signal for intermittency . This can be associated with a 
Renyi (multifractal) dimension 
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In the perturbative cascade these fluctuations in small bins are 
dominated by the possibility to have the tip of a jet within one 
bin. At asymptotic energies and large q we then get the result 
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□ partons, Qo=0.7GeV 
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FIG. 2: L3 results for ratios of factorial moments in cones around a 
jet axis, Fq{z)/Fq{Q), together with JETSET 7.4 PS predictions on 
partonic and hadronic levels. For details see ref. 0|. 
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FIG. 3: Analytic QCD predictions for A = 0.16 GeV. The curves 
show results for a.s = const., 3 different DLL approximations and a 
MLL result. For details see ref. tSJ. 



As for the distributions in the full phase space, there 
are large corrections from non-leading effects and from 
hadronization. Therefore the analytic results have had limited 
success in comparisons with experimental data. As examples 
figs. |2] and |3] show results from the L3 collaboration at LEP 
^ for factorial moments in regions separated by cones with 
opening angels ©o — © and ©o + © around a jet axis. The 
variable z is proportional to ln(©o/©). We see that although 
MC programs work quite well (fig. |2]i, the analytic calcula- 
tions do not reproduce the data (fig.©. We note that different 
NLL calculations give quite different results. These calcula- 
tions have different approximations for NNLL contributions, 
which shows that the higher order corrections cannot be ne- 
glected. (It was early realised that BE correlations give an 



important contribution to the multiplicity moments, and in the 
L3 analysis this effect is removed from the data.) 

It is also possible to define the A,-measure for a hadronic 
state instead of a parton state, and it appears that analytic re- 
sults for the A--measure are closer to the corresponding MC- 
results (j6|]. No comparisons with experimental data have been 
presented until now for such an analysis. 



D. Summary 

The gluon self coupling implies a fast growth of multiplic- 
ity and large fluctuations. In LLA we get the following asymp- 
totic results; 
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The QCD cascade has a fractal structure. The multiplicity 
in small phase space intervals shows intermittent features with 

a multifractal (Renyi) dimension ^ '2'Yq \ and the curve 

with the A,-measure has the dimension 1 +Yo ^. 

However, energy-momentum conservation and other non- 
leading effects are very important. NLL calculations give very 
different results if they differ at NNLL order. Hard emissions 
are the most important, but matrix elements factorize only for 
soft emissions. Therefore analytic calculations have in many 
cases met limited success. MC programs work, however, gen- 
erally quite well. 



II. FERMI-DIRAC CORRELATIONS IN pp AND AA PAIRS 
A. Momentum correlations 

The results on baryon coiTelations presented here are ob- 
tained in collaboration with R.M. Duran Delgado and L. 
Lonnbladd. 

Fermion correlation functions are usually fitted to the form 



CiQ)=N{l-Xe-'''Q'}. 



(18) 



The quantity R is here interpreted as the radius of the pro- 
duction region for the particle pair. The different LEP exper- 
iments have measured momentum correlations between AA 
and/or pp pairs, and obtained production radii around 0. 15fm, 
which is much smaller than the proton radius II ll Il2l Il3l |3 • 
In experiments the correlation function C{Q) is determined 
by comparing the observed number of pairs, N{Q), with a ref- 
erence sample, A^ref(2)- 



C{Q) 



m) 



(19) 



The important question is how to construct the reference sam- 
ple. Different methods have been applied: 
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• MC generation. This is model dependent. 

• Mixed events. The event structure is changed by gluon 
emission, which makes it difficult to construct a sample 
of events with similar structure. 

A common method to reduce the problems is to take a double 
ratio: 



C{Q) 



N{Q) I MC{Q) 



(20) 



Here the sample of Monte Carlo events, MC{Q), should be 
generated without including FD correlations in the model. It is 
argued that the effect of different event types should be similar 
in the real data and the MC events, and therefore be reduced 
in the ratio. 

The mixed pairs depend only on inclusive spectra, and the 
MC programs are tuned to reproduce the inclusive spectra 
with good precision. This implies that the tuning also gives 



which implies that 



NjQ) 
MC{Q) ■ 



(21) 



(22) 



Consequently the result depends very critically on a realistic 
MC, also if the double ratio is used in the analysis. 

In the Lund string hadronization model the colour string 
normally breaks by qq pair production. The ordering of the 
hadrons along the string, the rank ordering, agrees on average 
with the ordering in rapidity, with an average separation. Ay, 
of the order of half a unit in rapidity. A baryon-antibaryon pair 
can be produced when a breakup is generated by a diquark- 
antidiquark pair, with quark-antiquark pairs on either side. In 
this picture two baryons must always be separated by at least 
one antibaryon (and normally also with one or more mesons). 
This will give a strong anti-correlation between two baryons 
in rapidity and in momentum. This is seen in fig. |4] which 
shows correlations in pp and pp pairs. We see that there is 
a positive correlation between protons and antiprotons, which 
are neighbours in rank, but a negative correlation between two 
protons. 

In the model there is also a correlation between momen- 
tum and space coordinates for the produced hadrons. Thus 
two identical baryons are (in the model) well separated also 
in coordinate space, and we would from this picture expect 
Fermi-Dirac correlations to correspond to a radius ~ 2 — 3 fm. 

We see that the range for the pp correlation in fig.|4]is given 
by 2 ^ 1.5GeV ^1/ (O.lSfm). We note that this corresponds 
exactly to the correlation length reported in the experiments, 
although in this model the production radius is very much 
larger. The strength of the correlation is, however, smaller 
in the MC than in the data. 

This raises the question: Is the difference between data and 
MC really a FD effect, or could the MC underestimate the 
strength of the correlation? 

There are a number of sources for uncertainty in the MC: 
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FIG. 4: MC results for the ratio Cmc(2) MC{Q)/MCmix{Q) for 
pp- (x) and pp-paus (+). 



MC proton correlations 




a=0.1 b=0.45 GeV" 

a=0.3 b=0.58 GeV" 

a=0.354 b=0.523 GeV" 

a-0.5 b-0.7 GeV 

a-1 b-1 GeV" 

a^2 , b^2 , GeV 

1 7 8 



Q (GeV) 



FIG. 5: The ratio Cmc{Q) = MC{Q)/MCmix{Q) for different val- 
ues of the parameters a and b. Larger (a,Z7)-values give a stronger 
correlation and a deeper dip for small Q. 



1) There are two fundamental parameters, a and b, in the 
Lund model "splitting function" 



(23) 



The hadron multiplicity depends essentially on the ratio (a + 
l)/b. This ratio is therefore well determined by experiments, 
but a and b separately are more uncertain. Small values of 
a and b correspond to a wide distribution f{z), and a wide 
distribution in the separation. Ay, between hadrons which are 
neighbours in rank. Large values of a and b imply a narrow 
Ay- distribution and therefore lower probability for two parti- 
cles to be close in momentum space. The effect of varying a 
and b keeping the multiplicity unchanged is shown in fig.|5] 

2) The parameter bis a universal constant, but a may be dif- 
ferent for baryons, although data are well fitted by a universal 
fl-value. This gives some extra uncertainty. 

3) The MC does not exactly reproduce the Lund hadroniza- 
tion model. The splitting function in eq. $23[ gives a correct 
result when the remaining energy in the system is large. To 
minimise the error at the end of the cascade, when the energy 
is small, the MC cuts off hadrons from both ends randomly, 
and joins the two ends when the remaining mass is small. This 
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FIG. 6: A single jet without a junction has a large dip in C^ciQ) = 
MC{Q) /MCjnixiQ) fo'' small 2 (x). This correlation is reduced by 
the approximate treatment in the MC of the small mass systems close 
to the "junction" (+). 
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FIG. 7: Cmc{Q) = MC{Q) /MC^i^iQ) with (+) and without (x) a 
15% admixture of pions. 



works well for inclusive distributions, but implies that the cor- 
relations do not correspond to the model prediction for parti- 
cles close to the "junction". Fig.|6]shows results for a single 
jet without a junction compared with the standard result. 

4) Gluon emissions imply that straight string pieces are 
small compared to the mass of a BBB system. This gives also 



extra uncertainty. 

5) Pion correlations are not perfectly reproduced by the 
MC. As an example the DELPHI antiproton sample contains 
15% pions. Fig. [T] shows results with and without a 15% 
pion admixture. We see that an error in the simulation of 
the pion-pion or pion-proton correlations also affects the esti- 
mated proton-proton correlations. 

In summary we see that there are many effects which make 
the MC predictions for pp or AA correlations quite uncertain. 
As the experimental determination of the correlations rely so 
strongly on a correct MC, it is therfore at present premature to 
conclude that the production radius has the very small value 
around 0.1 5fm. 

B. Spin-spin correlations 

A particles reveal their spin in the orientation of their decay 
products. This has been used to study AA correlations without 
the need for a comparison with MC results. A AA pair with 
total spin 1 must have an antisymmetric spacial wave function 
and is therefore expected to show a suppression for small rel- 
ative momenta Q. AA pairs with total spin has a symmetric 
spacial wave function, and should therefore show an enhance- 
ment for small Q, similar to the correlation for bosons. There- 
fore one expects pairs with small Q-values to be dominantly 
5 = 0. Analyses at LEP ful [H Il6f l do indicate such an ef- 
fect. They are however based on rather low statistics, and the 
errors are presently too large for any definite conclusion about 
the strength and range of the effect. 



C. Summary 

A production radius R ^ 0. 15fm for baryon pairs is not con- 
sistent with the conventional picture of string fragmentation. 

Experimental results for pp and AA correlations depend 
sensitively on a reliable MC. 

The observed (anti-)correlation has the same range in Q as 
the correlation in the MC, but is stronger 

Uncertainties in the MC implementation are large. 

Conclusion: It is therefore premature to claim evidence for 
a new production mechanism. 
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